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Abstract 

We obtain large deviation results for a two time-scale model of jump-diffusion processes. 
The processes on the two time scales are fully inter-dependent, the slow process has small 
perturbative noise and the fast process is ergodic. Our results extend previous large devi¬ 
ation results for diffusions. We provide concrete examples in their applications to finance 
and biology, with an explicit calculation of the large deviation rate function. 
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1 Introduction 


For a number of processes in finance and biology the appropriate stochastic modelling is done in 
terms of multi-scale Markov processes with fully dependent slow and fast fluctuating variables. 
The most common examples of such multi-scale processes (random evolutions, diffusions, state 
dependent Markov chains) are all particular cases of jump-diffusions. The law of large numbers 
limit, central limit theorem, and the corresponding large deviations behaviour of these models 
are all of interest in applications. 

One case of their use in finance is in multi-factor stochastic volatility models, which are 
used to capture the smiles and skews of implied volatility. The separation of time scales is 
helpful for calibration, since it allows one to reduce the number of group parameters. The rate 
function from the large deviation principle for the stock price process can be used to obtain 
the price of short maturity options, as well as the limit of the at-the-money implied volatility. 
These have been explicitly calculated for models in which the logarithm of the stock price and 
the stochastic volatility are driven by diffusions 1 [FFK12] . [FFFIO] ). However, much of the 
empirical evidence f |B-NS01] . |Kou02] 1 suggests that mean-reverting jump-diffusions would be 
a more appropriate model for the problem. 

In biology one case of their use is in models of intracellular biochemical reactions. Due to 
low copy numbers of various key molecular types and varying strengths in chemical bonds, nor¬ 
malized copy numbers of different types of molecules are processes on multiple time-scales (see 
|BKRP05] , [KK13] for references to the biology literature). Changes in molecular compositions 
are modelled by state-dependent Markov chains, and on the slower time scale are well approxi¬ 
mated by diffusions with small noise or piecewise deterministic Markov chains ( [KKP12| ). The 
rate function from the large deviation principle for slowly fluctuating molecular species is used 
to calculate the propensity for switching in a network that has multiple stable equilibria. Since 
intracellular processes are also subject to other sources of ‘extrinsic’ noise, multiple time-scale 
diffusions may include jumps from additional sources. For example, there can be errors during 
cell division (!P au m, [p au m-, a stochastic model combining both reactions and cell division 
was analyzed in [McSP14] . 

Large deviation results for multi-scale diffusions have been studied by Freidlin (see [F98j 
Chapter 7), Veretennikov ( [VerOO] ). Dupuis et al. ( |Dupl2| ), and Puhalskii f |Puhl5| l. For the 
multi-scale Markov chains where the slow process is a piecewise deterministic Markov processes 
and the fast process is a Markov chain on a finite state space explicit results were obtained by 
Faggionato et al. ( |Fagg09 FagglOj ). For jump-diffusions there are very few large deviation 
results. On a single time scale, there are results by Imkeller et al. ( [Imkn9] ) for first exit times 
for SDEs driven by symmetric stable and exponentially light-tail symmetric Levy processes. 
An approach based on control theory and the variational representation was developed by 
Budhiraja et al. in [Budll] and extended to infinite dimensional versions |Budl3] (that is, 
SPDEs rather than SDEs driven by a Poisson random measure). It is not easy to see how 
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to use these results in a multi-scale model of jump-diffusions. A special case of a multi-scale 
process where the slow process is a diffusion and the fast process is a mean-reverting process 
driven by a Levy process was studied by Bardi et al. f |BCSj L and the authors use PDE 
methods to prove asymptotics of an optimal control problem. 

A general method for Markov processes based on non-linear semigroups and viscosity meth¬ 
ods was developed by Feng and Kurtz in [FK06] . However, verifying the abstract conditions 
needed to apply this method to multi-scale jump-diffusions is a non-trivial task. In this pa¬ 
per we give a proof of large deviations for two time-scale jump-diffusions, using a technique 
developed by Feng et al. in |FFK12| . The advantage of this method is that it is constructive 
and, with some effort, can be tailored to different multi-scale processes. Our proof follows the 
steps of |FFK12| , extending it to processes with jumps and full dependence of the slow and fast 
components. It is based on viscosity solutions to the Cauchy problem for a sequence of partial 
integro-differential equations and uses a construction of the sub- and super-solutions to related 
Cauchy problems as in [FFK12j . Our results hold for slow and fast jump-diffusions which are 
fully inter-dependent, and where the fast processes is ergodic but not necessarily symmetric. 
In case the evolution of both processes is spatially homogeneous in the slow variables, we can 
also provide a more explicit (than a solution to a variational problem) formula for the rate 
function. 

2 Two time-scale jump-diffusion 

Consider the following system of stochastic differential equations: 

dx, ,t =b{X,,t-,ye,t-)dt + ebo{X,,t-,ye,t-)dt + 

+ €j k{X^^t-,ye,t-,z)N^'^^\dz,dt), 

dy, ,t =\{x,,t-,ye,t-)dt + 

+ J h{x,^t-,ye,t-,z)m-(^\dz,dt), 

y^e,0 = 3 ^ 0 , ye,0 = UO, 

where iVe ■), N, ■) are independent Poisson random measures with intensity measures 
vi{dz) X ^dt,V 2 {dz) X ^dt; the Levy measures ui and 1^2 satisfy Jg(I A z‘^)v 2 {dz) < 00 and 
/jg(l A z^)v 2 {dz) < 00 ; the centered versions are defined as 

_/V«•) = •) — ui{dz) X -dt, •) = •) — V2{dz) x -dt 

and are independent Brownian motions independent of N 

To ensure existence and uniqueness of solutions to the system m we assume 


(la) 

(lb) 
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Assumption 2.1 (Lipschitz condition). There exists Ki> Q such that V(xi,yi), (x 2 ,?/ 2 ) £ 

\b{x2,y2) - h{xi,yi)\^ + \ho{x2,y2) - ho{xi,yi)\^ + \hi{x2,y2) - 6i(xi,yi)p 
+ \<y{x2,y2) - o-(xi,yi)p + |<Ti(x2,?/2) “ CFi{xi,yi)\’^ 

|fc(x2,y2,^;) - k{xi,yi,z)\^vi{z)dz + / \ki{x2,y2,z) - ki{xi,yi,z)\^U2{z)dz 

< Ki{\x2 - Xil^ + \y2 - yi\^). 

Assumption 2.2 (Growth condition). There exists K 2 > 0 such that\/{x,y) E 
\b{x,y)\‘^ + \bo{x,y)\‘^ + |6i(x,7/)p + \cr{x,y)\‘^ + \ai{x,y)\‘^ 


+ / \ki{x,y,z)\‘^U 2 {z)dz + / \k{x,y, z)\‘^iyi{z)dz < K 2{1 + x'^ + y"^). 


Define 


V{y; x,p) := b{x, y)p + y)p^ + J - 1 - pk{x, y, z)^ ui{z)dz. 

For each x and p in U. there exists Kx,p > —00 such that 

V{y]x,p) > Vy E M. 


(3) 


(4) 


(5) 


If existence and uniqueness of solutions to (fTa1) + (llbp can be established by other means, we 
will only assume the growth condition i.e. Assumption 12.21 that the coefficients are continuous, 
and the lower bound ([5]) on V. 

The infinitesimal generator of (Ag,!^) is for / E C'|(M x M) defined by 
k^ef(x, y) = b{x, y)da;f{x, y) + pa{x, y)ai (x, y)dlyf{x, y) 


+ ebo{x,y)dccf{x,y) + -a^{x,y)dl^f{x,y) 

+ - I ifix + ek{x,y,z),y)-f{x,y)-ek{x,y,z)dxf{x,y))ui{z)dz 


+ 


( 6 ) 


bi{x,y)dyf{x,y) + -al{x,y)dyyf{x,y) 


+ J {f{x,y + ki{x,y,z)) - f{x,y) - ki{x,y,z)dyf{x,y))u 2 iz)dz 
Fix X E M and let denote the process satisfying the SDE 
dYt =bi{x,Yt-)dt + ai(x,yt_) 

+ [ ki{x,Yt-,z)N^‘^\dz,dt), Yq = yo. 


(7) 


This is the SDE (flbl) where e is set equal to 1 and is set equal to x. Let Cf denote the 
generator of Y^, then, for / E C^M), 


71f/(y) :=bi{x,y)dyf{y) + -a‘l{x,y)dlyf{y) 

+ / {f{y + ki{x,y,z))-f{y)-ki{x,y,z)dyf{x,y))u 2 {z)dz. 


( 8 ) 
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For fixed p E M define the perturbed generator for / E by 


: = [pcr{x,y)(Ti{x,y)p + bi{x,y)]dyf{y) + -al{x,y)dyyf{y) 

+ / if{y + kiix,y,z))-f{y)-ki{x,y,z)dyf{x,y))u 2 {z)dz, 


( 9 ) 


and let Y^’P be the process corresponding to the generator For each E M we assume 
the following about Y^’P 

Assumption 2.3 (Ergodicity condition). The process Y^’P is Feller continuous with transition 
probability p^’^{yo,dy), which at t = 1 has a positive density Pi’^{yo,y) with respect to some 
reference measure a{dy). 

Assumption 2.4 (Lyapunov condition). There exists a positive function ('(•) E C^(M), such 
that ( has compact finite level sets, and for each compact set F C K, 0 E (0,1] and I E M, there 
exists a compact set C M such that 

{y E M : -ee~‘^£'l'^e‘^{y) - {\V{y,x,p)\ + \bo{x,y)p\ +a‘^{x,y)) < 1} C Ai^g^r, Vp E F,Vx E M. 

( 10 ) 


Remark 2.1. In the case where the domain ofY is compact, we can define C = 0 which will 
satisfy Assumption \2.4\ 


Remark 2.2. Some arguments are simpler in the special case Y^’P in addition has a unique 
invariant probability measure ttP{x, •) with respect to which pf’^(yo, y) is symmetric and 7rP{x, ■) 
is reversible, that is 


k^i’^fiyWix,y)dy 


0, V/ E C~(M). 


and 

j f{y)k^l’^ 9 {y)T^^{x,y)dy = j g{y)Cl''^f{y)T:P{x,y)dy, V/,yEC'^(M) 


2.1 Examples 

We give some examples of Y that satisfy Assumption 12.31 as well as a multiplicative ergodicity 
condition of the form 

e~'^Cl''^e^{y) < -C,{y) + d 

for Q with compact level sets and some constant d > 0. One needs to know the coefficients of 
the process X to know whether these examples also satisfy Assumption 12.41 Define VP{x, y) := 
V{y;x,p) + |6o(x,y)| + a‘^{x,y). If VP{-,-) is a bounded function for bounded p, then the 
multiplicative ergodicity condition is sufficient for Assumption 2.4 to hold. If VP{x,y) is an 
unbounded function but has compact level sets, and if the E-multiplicative ergodicity condition 
of the form 

e“‘’£^’^e^(y) < —cVP{x, y) + d, for some c > 1, d > 0 

is met for C, with compact hnite level sets, then it may be possible to use this condition in place 
of Assumption 12.41 and obtain all the same results (see Example 14.11 and Remark 14.Ij) . 
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Example 2.1. Let p = 0, bi(x, y) = —hi(x)y, cti(x, y) = cri{x) and ki(x, y, z) = — y, 

y/bi{x} 

where bi{x), ai{x) > 0 are continuous. Let 1 ^ 2 ( 2 ) = exp {—Since the intensity measure ui 
is a bounded measure, we use instead of the compensated Poisson process . For each 
X G M, the solution to 

dYf = -bi{x)Yfdt + (Ti(x)dlEf ^ + [ (^^kLz - Yf)N^^\dz, dt) 

Jr-{o} ■\/bi{x) 


has unique invariant probability distribution 7r{x, dy) = y^|^exp{— 
symmetric with respect to it. Geometric ergodicity is satisfied by C{y) := 


bi{x)y'^ 

Ffix) 

bi{x) 

2F[(x) 



and Y^ is 


Example 2.2. Take p = 0 and let a G (1)2). Let Zt be a 1-dimensional symmetric Levy 
process whose Levy measure is V 2 {z)dz = \z\~^^~^°^'>l\z\yidz. Its infinitesimal generator is the 

truncated fractional Laplacian —(—defined as 


-i-ApTm 


j (/fe + z) - 


for /gC2(R). 


Let ai{x,y) := a{x)ai{y) where a(-),cJi(-) > 0 are such that a(-) is continuous and fTi(-) is 
locally l/a-Holder continuous and liminf|y|_).cc > 0- 


dYfi = afix,Yf_)dZt. 


Then from Theorem 1.7(i) in fCWlfij , Tr{x, dy) := unique invariant probability 

measure for the Y^ process and Y^ is tt{x, )-reversible. From Lemma 3.2 in icwm , we get 
C{y) := ln(l + \y\^) for 6 G (0,1) satisfies the geometric ergodicity condition. The special case 
of this example with ai = l is also considered in 

Example 2.3. Let c(z,z') be a non-symmetric function such that 0 < cq < c(z,z') < ci, 
c(z, z') = c(z, —z') and \c{z, z") — c{z', z")\ < C 2 I 2 : — z'jl^ for some fi G (0,1). Let a G (0, 2), 
and Zt be a 1-dimensional non-symmetric process whose infinitesimal generator is defined by 

k^cfiy) = [ ifiy + ^) - fiy) f^^i'n+a^^ dz, for f G 


Let 


dYt = -Ytdt + dZt. 


Heat kernel estimates from \CZ13^ imply this non-symmetric jump diffusion is Feller continuous 
with a positive transition density pt{yo,y),\/t > 0 


Example 2.4. Let Y^ be a birth-death Markov chain with birth rate r^{y) = A(x) and death 
rate r^fy) = p{x)y, satisfying \{x),p{x) > 0. Since its state space is countable its transition 
density is positive, with a unique reversible invariant distribution tt{x, y) = ^ 

2/G{0,1,...}. 
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3 Large deviation principle 


We prove a large deviation principle for as e —0 using the viscosity solution approach 

to verify convergence of a sequence of exponential generators. Define 


u^{t,x,y) := elnE 


e ^ \X^^Q = x,Y^^Q = y 


( 11 ) 


where h E C'ft(M), the space of bounded uniformly continuous functions on M. It can be shown 
(see [FK06] ) that for each h E C'b(M), is a viscosity solution of the Cauchy problem: 


dtu = H^u in (0, T] x M x M, 
u(0, X, y) = h{x), for (x, y) E K x M, 

where the non-linear operator is the exponential generator: 

H,u{x,y) := 

= b{x, y)da:u{x, y) + pa{x, y)ai{x, y)dlyu{x, y) + y){da;u{x, y)f 


( 12 ) 


+ e 


bo{x,y)dj:u{x,y) + l-a‘^{x,y)d^^u{x,y) 


u(x + ek(x,v,z),v)—u(x,v) . 

e J - 1 - k{x, y, z)dxu{x, y) vi{z)dz 


1 

+ - 
e L' 




1 


pcr(x, y)ai{x, y)dxu{x, y)dyu{x, y) 6i(x, y)dyu{x, y) + ^^a^x, y)dyyu{x, y) 


e ^ - 1 - dyu{x, y) \ i^ 2 {z)dz + -^aj{x,y){dyu{x,y)f. 


(13) 


In systems with averaging under the law of large number scaling we can identify the limiting 
non-linear operator Hq as the solution to an eigenvalue problem for the driving process 
obtained from Yf, with = x and e = 1. 

We first identify uq, the limit of Ue as e —>■ 0, using heuristic arguments. Assume 

Ue{t, X, y) = uo{t, x) + eui{t, x, y) + e^U 2 {t, x,y) + .... (14) 

Using the e expansion of m, in equation (fT^ . and collecting terms of 0(1), we get 

dtUo{t,x) = b{x,y)dxUo{t,x) + ^a‘^{x,y){dxUo{t,x)f 

+ J t,i{z)dz 

+ P<x{x,y)(Ji{x,y)da;Uo{t,x)dyUi{t,x,y) -h bi{x,y)dyUi{t,x,y) + ^al{x,y)dyyUi{t,x,y) 
+ i^(^i{x,y){dyUi{t,x,y)f. 

(15) 
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Please note that as this is merely a formal derivation, we have ignored some technical details 
(snch as justifying interchanging the limit and integral to get the the second line in the above 
equation). The rigorous proof that follows shows that this formal derivation is indeed correct. 
Denote dxUo{t,x) by p and dtUo{t,x) by A. Fix t,x and hence p and A. Using the perturbed 
Cl generator ([U]), the equation (fTCD can be written as an eigenvalue problem: 

(C'l'P+ V{y;x,p))e^^ = \e^\ (16) 

where V is as defined in O- Note that the eigenvalue A depends on x and p, and that if we 
write Hq{x,p) := A then uq satisfies 

dtUo{t,x) = Ho{x,dxUoit,x)), 

In the rigorous proof that follows, we identify the limiting operator Hq to he as defined in (jl8h 
which is shown in |DV75] to be the principal eigenvalue A in (|16p . By the expansion (I14p . it is 
clear that tto(0, x) = h{x). 

The approach of [FK06| for obtaining the large deviation principle is to prove convergence 
of nonlinear semigroups associated with the nonlinear operators Fig. In [FK06| the first step 
is identifying the limit operator Hq. Existence and uniqueness of the limiting semigroup is 
obtained by verifying the ‘range condition’ for the limit operator. This amounts to showing 
existence of solutions to the equation (I — aHo)/ = h for small enough a > 0 and sufficiently 
large class of functions h. Since the range condition is difficult to verify, a viscosity method 
approach is adopted and the range condition is replaced with a comparison principle condi¬ 
tion for (I — aHo)f = h. In the viscosity method, existence of the limiting semigroup is by 
construction, while uniqueness is obtained via the comparison principle. 

The approach in this paper uses convergence of viscosity solutions to the Cauchy problem 
for PIDEs (II2p . and to show existence and uniqueness of the limit one then needs to verify the 
comparison principle for the Cauchy problem dtUo{t,x) = Ho{x, dxUo{t, x)), with uo{0,x) = 
h{x). 

In the proof of the comparison principle we will also use a Donsker-Varadhan variational rep¬ 
resentation l [DV75p for Hq as follows. Let P(M) denote the space of probability measures on 
M. Define the rate function J{-]x,p) : P(M) i-)- M U {-|-oo} by 

J{p;x,p):=— inf / —(17) 

geD++{c^p^) 9 

where (C^’^) C C'b(M) denotes the domain of C^'^ with functions that are strictly bounded 
below by a positive constant. Then |DV75j implies that the principal eigenvalue Ho{x,p) = A 
in (fT6]l is also given by 

Ho{x,p)= sup { V{y;x,p)dp{y)-J{p-x,p)], (18) 

\J J 

where V{y;x,p) = b{x,y)p + \a‘^{x,y)p^ + f — l) vi{z)dz. 













Remark 3.1. In the special case also has a reversible invariant measure 7r^(x,), we can 
use the Dirichlet form representation for J. Define the Dirichlet form associated with by 


Then, Theorem l.ff in Stroock fStrooBf}/ implies that 


J{fj.]x,p) = 


(V \ldJ^) < 7rP(x, •) 

if hi-) 


+00 


The variational formula il8\) then reduces to the classical Rayleigh-Ritz formula 

Ho{x,p)= sup ([viy-,x,p)fiy)d7rP{x,y)dy + {C'l’Pf,f)\. 

/GL2(7rP),|/P = l \J J 


(19) 


( 20 ) 


To sum up, we will prove that: 

Lemma 1. Let Hq be as defined in dm), and suppose the comparison principle holds for the 
nonlinear Cauchy problem: 

dtUo{t,x) = Ho{x,dxUo{t,x)), for {t,x) G (0,r] x M; 

Mo(0, x) = /l(x). 

Under the Assumvtions \2.1^2.4\ the sequence of functions {u^}e>o defined in m converges 
uniformly over compact subsets of [0,T] x M x M as e —>■ 0 to the unique continuous viscosity 
solution Uq of m- 

Lemma 2. The sequence of processes {Xe^i}e>o is exponentially tight. 

Theorem 3. Let suppose all the Assumptions from Lemma [I] hold. Then, 

{Xf:^t}e>o satisfies a large deviation principle with speed 1/e and good rate function 

I{x,xo,t)= sup {/i(x) -'Uo(Lxo)}. (22) 

/iGCb(R) 

Proof. By Bryc’s theorem (Theorem 4.4.2 in [DeZ98] l. Lemmas [T] and [2] give us a large deviation 
principle for {Xe_t}e>o as e —0 with speed 1/e and good rate function I given by ([2^ . □ 

One of the key conditions for Lemma [T] requires one to check that the comparison principle 
holds for Hq. This condition cannot be established using only the general Assumptions 12.11 
[M and needs to be verified on a case by case basis. However, standard theory of comparison 
principles for viscosity solutions (Theorem 3.7 and Remark 3.8 in Chapter II of |BD97j ) implies 
that it does hold for m as soon as Hq is uniformly continuous in x,p on compact sets (see 
Lemma fTOl of the Appendix). In some cases Hq can be explicitly calculated (see Example 14. 2 p 
and continuity directly verified. In other cases one may need to resort to proving that the 
expression as on the right-hand side of (IA.37P is non-positive, using the specifics for the case 
at hand. 
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Corollary 4. Any of the following separate sets of conditions are sufficient for the comparison 
principle for the non-linear Cauchy problem [2J\) to hold: 

(i) Hq is uniformly continuous in x,p on compact sets; 

(a) the eoeffieientsbi{x,y),ai{x,y),ki{x,y,z) are independent of x, the coefficients b{x,-), a(x,-) 
are bounded (bounded functions of y) for each x, and p = 0 i.e. the correlation between 
the Brownian motions driving X and Y is 0. 

Proof. For (i) see Lemma[iniof the Appendix which is based on Theorem 3.7 and Remark 3.8 
in Chapter II of |BD97] . 

For (a) we can directly verify that under these conditions Hq{x,p), given in (I18p . is uniformly 
continuous on compact sets of x and p. For this, first observe that under the conditions in 
(a) the rate function J in (fT8]l will be independent of x and p. Additionally, f V{y, x,p)dfr{y) 
is uniformly Lipschitz in x and p (uniform over all p, G V(M.) ), over compact sets of x and 
p. Finally, since the supremum of uniformly Lipschitz functions is uniformly continuous over 
compact sets, we have the result. □ 

Note that in Corollary |U condition (i) is a more general condition and (ii) is a sufficient 
condition (on the coefficients of the model) under which condition (i) holds. 

In very special cases, we can further simplify the expression for the rate function: 

Corollary 5. If the coefficients in the SDE ([I]) are independent of x, then Hq{x,p) becomes 
Hq{p) and by Lemma D.l in ]FFKlS!f . we get 

I{x,-xo,t) = tLo ^ , (23) 


where Lo(') is the Legendre transform of Hof). 


The proof of Lemma [U takes up the bulk of the paper, and consists of the following steps. 

(Sec 13.ip • By taking appropriate limits of solutions to the Cauchy problem (fl^ we construct 
upper-semicontinuous and lower-semicontinuous functions and yf, respectively; 

• Using an indexing set a G A, we construct a family of operators Hof ; a) and Hif ; a), 
in such a way that the upper-semicontinuous function is a subsolution to the Cauchy 
problem for the operator infagA{77o(' i®)}) and the lower-semicontinuous function yf is 
a supersolution to the Cauchy problem for the operator sup^^/^{Hif ; a)}. 

fSec 13.2p • We prove a comparison principle between subsolutions of infQ,gA{77o('; «)} and super¬ 
solutions sup„gA{^i(' ^'a)} above; 

• We show that this comparison principle implies convergence of solutions uf to the 
Cauchy problem (fT^ for H^ to solutions Uq to the Cauchy problem (f^ for Ho. 

The proof of Lemma [2] uses the estimates obtained in the proof of Lemma [U (Section 13.3|) . 
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3.1 Convergence of viscosity solutions of PIDEs 

In Lemma [T] we use notions of viscosity solutions, subsolutions and supersolutions. For the 
standard meaning of these terms, as well as for the definition of the comparison principle, we 
refer the reader to Definition 4.1 in |FFK12j . Their extension to partial integro-differential 
equations (PIDEs) was obtained already in [Alv96| and can be found in [Bar08] . 

The proof of convergence of to Uq follows the same steps as Lemma 4.1 in |FFK12j which 
carries over directly to viscosity solutions of PIDEs. Because we will need to verify that the 
conditions there are met, we restate Lemma 4.1 from |EEK12] for viscosity solutions of PIDEs. 

Let {fff;}£>o denote a family of integro-differential operators defined on the domain of 
functions l)+ U Z)_ where 

D+ := {/ : / G lim inf f(z) = -hoo} 

r->-oo |2|>r 

D- := { —/ : / G lim inf f{z) = +oo}. 

r^oo \z\>r 

Define domains analogously replacing by M. Consider a class of compact sets in 

M X M defined by 

Q := {K X K : compact K,K CC M}. 

Let Ug be the viscosity solution of the Cauchy problem dtu = H^u for the above operator LIg, 
with initial value h, and define 

Definition 3.1. 

ui{t,x) := sup{limsupUg (L,Xg,?/g) : 3(L,Xg,7/g) G [0,r] x iL x K, 

e^0+ 

(L,Xg) ^ it,x),K X K G Q}, 

Ui{t, x) := inf{liminf (L, Xg, y^) : 3(tg, Xg, j/g) G [0, T] x AT x K, 

(tg,Xg) ^ it,x),K X K G Q}. 

Define to be the upper semicontinuous regularization ofup andyfi the lower semicontinuous 
regularization ofu^. 

Einally, define the limiting operators (which will be first-order differential operators) Hq, Hi 
on domains D_|_ and Zl_ respectively, as follows. Let A be some indexing set, and 

Hi{x,p] a) : M X M !-)■ M, a G A, z = 0,1. 

Define Hof{x) := HQ{x,dxf{x)), for / G and Hif{x) := Hi{x,dxf{x), for / G Zl_, where 

Ho{x,p) := inf Ho{x,p]a), 
qGA 

Hi{x,p) := sup Hi{x,p;a). 
qGA 

Henceforth, with slight abuse of notation, we will refer to Hfi-, •) as operators. 

Suppose the following conditions hold: 
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Condition 3.1 (limsup convergence of operators). For each /o G -D+ and a ^ A, there exists 
/o,e G (which may depend on a) such that 

1. for each c > 0, there exists K x K ^ Q satisfying 

{{x,y) : > -c} n {{x,y) : /o,e(x,?/) < c} C K x K] 

2. for each K x K G Q, 

lim sup |/o,e(x,y) -/o(x)| = 0; (24) 

{x,y)£KxK 

3. whenever {x^^yf) ^ K x K ^ Q satisfies x^ x, 

limsupFe/o,e(xe,ye) < V/o(x); a). (25) 

£->■0 

Condition 3.2 (liminf convergence of operators). For each fi G D- and a G A, there exists 
/i^e G -D_ (which may depend on a) such that 

1. for each c > 0, there exists K x K ^ Q satisfying 

{{x,y) : HJi^^{x,y) < c} n {{x,y) : fi,eix,y) > -c} C K x K] 

2. for each K x K ^ Q, 

lim sup |/i(x) -/i,e(x,y)| = 0; 

{x,y)eKxK 

3. whenever {xe,ye) £ K x K G Q, and Xe x, 

Wvami HJi^^{x^,y^) > Hi{x,Vfi{x);a). 

e^O 

In this case the following convergence results for as e —^ 0 hold. 

Lemma 6. Suppose the viscosity solutions to the partial integro-differential equation 

dtu = H^u, u(0,x) = h(x) 

are uniformly bounded, sup^>Q ||u(ll| < oo. Then, under Condition \3.1{ uf is a subsolution of 

dtu{t,x) < HQ{x,'\/u{t,x)) (26) 

and, under Condition \3.2i yf is a supersolution of 

dtu{t,x) > Hi{x,'Vu{t,x)). (27) 

with the same initial conditions. 
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As the proof is the same as the proof of Lemma 4.1 in [FFK12| we omit it here. We do need to 
check Conditions 13.11 and 13.21 hold for our problem. This involves identifying the right indexing 
set A, the family of operators and and the appropriate test functions /o,e 

and /i^£, for each given /o and /i, respectively. 


Verifying Condition I3.lt As in [FFK12] . we let 

A:={(e,0):CGC2(R),O<0<l} 


and we define the sequence of operators as in (1131) on the domain 

D+ ■= if € : f{x) = 0(x) +7log(l + G C'c^(M),7 > 0}. 

Define the family of operators Hq {x,p]^,e) for (^,6») G A by 

Hoix,p-, 9) := sup{6(x, y)p + y)p^ + f _ i _ z)) vi{z)dz 

y&M. Z J \ / 

+ (1 - 6)e~^C\'^e^{y) + 6e~^Cl''^e^{y)}. 

For any / G Dj^ and (^, 0) G A define a sequence of functions 


(28) 


h,e{x,y) ■.= f{x) + eg{y), where g{y) ■.= {I - e)i{y) + OQiy), 
and C is the Lyapunov function on M satisfying Assumption 12.41 Then, 


1 


Hefo,e{x, y) = b{x, y)d^/(x) + -a (x, y){d^f{x)) + e bo{x, y)d^fix) + -a (x, y)d^^f{x 


1 


+ 


f {x-\-€k{x ,y ,z) ,y)— f {x ,y) 


1 


- 1 - k{x,y,z)dxf{x)^ ui{z)dz + e 


< b{x, y)d^f{x) + -a {x, y){d^f{x)) + e 6 o(x, y)d^f{x) + -a {x, y)d^^f{x 


1 


+ 


f (x-\-ek(x .y ,z^ .y^ — f(x.y) . 

e s - 1 - k{x, y, z)dj,f{x) vi{z)dz 


+ (1 - + ee-^cf-f^^'^e^{y). 

so, for any sequence {xe,ye) such that x^^x 


(29) 


\misxvp HjQ^^{x^,y^) < Ho{x,d^f{x)-,C,d), 

e^O 

thus verifying Condition 13.1131 holds. 

By choice of D+, / G has compact level sets in R. Also note that ||dx/|| + ||5*a;/|| < oo. 
Assumption [23] ensures that — FFe/o,e(x, •) has compact level sets for all x in compact sets. This 
proves Condition 13. 1111 holds. Condition 13.1121 is obvious by choice of functions /o,^. 


Verifying Condition 13.21 is exactly the same as verifying Condition 13.11 except that the 
sequence of operators FFg are now defined on the domain 

D_ := {/ G ^^(R) : f{x) = (p{x) - 7 log(l + x^); (j) G > 0}; 
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the family of operators Hi{x,p]^,6) for {^,0) G A is defined by 

Hi{x,p-, 6) := inf {6(x, y)p + y)p^ + j - 1 - pk{x, y, z)^ vi{z)dz 

+ (1 + e)e~^Cl’^e^{y) - ee~‘^Cl'^e^{y)}] 

and for any / G Z?_ and 0 G A the sequence is defined as 

fiA^, y) ■= fix) + €g{y), for g{y) := (1 + 6»)^(y) - 6»C(?/), 
so that for any sequence (x^, y^) such that x^ ^ x we now have 


limini HJiAxe,ye) > Hiix,d^f{x);^,e) 

£->■0 

verifies Condition l.S, 2131 holds. Conditions 13.2111 and l3.2l2l hold by the same arguments as above. 

3.2 Comparison Principle 

The rest of the claim of Lemma [T] requires proving uniqueness of solutions to dtu = Hqu, with 
initial value h. This can be verified using the comparison principle on the subsolutions and 
supersolutions of the constructed limiting operators Ho and Hi, and the variational represen¬ 
tation of Ho from (fT8|) . We use the following Lemma 4.2 from |FFK12| . 

Lemma 7. Let yf and lA he defined as in Definition \3.1{ If a comparison principle between 
subsolutions of (fM]l and supersolutions of (f27l) holds, that is, if every subsolution vi of (f26l) 
and every supersolution V 2 of (EZD satisfy vi < V 2 , then yfi = and uAt,x,y) — UQ{t,x), 
where Uq := yfi = lA, as e ^ 0, uniformly over compact subsets of [0,T] x M x R. 

Proof. The comparison principle gives ^ vf, while by construction we have ufi <tA. This 
gives uniform convergence of —)• uq := = lA over compact subsets of [0, T] x R x R. □ 

We next prove the comparison principle for subsolutions of (|M]l and supersolutions of (|27|l . 
that is every sub solution of 

dtu{t,x) < Ho{x,p) ■.= inf Hoix,p;^,9), 

0<e<i,5eC2(K) 

where Hq is as defined in (I28p . is less than or equal to every super solution of 

dtuit,x) > Hiix,p) := sup Hi{x,p;f,,9) 

o<e<i,5ec2(R) 

where Hi is as defined in (1301) . We follow the steps in Section 5.2 in |FFK12j with some 
modifications. The key step is proving 

Operator Inequality: 

inf Hoix,p-,9,f) < Hoix,p) < sup Hi{x,p-,9,Q, (31) 

o<0<i,§ec2(R) o<e<i,5ec2(R) 

where Ho{x,p) is as defined in (fTHI) . 
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Recall the definition of the rate function J from (I17p and variational representation ol Hq as 


Ho{x,p) = sup 
mg'P(R) 


V{y;x,p)dn{y) - J{y,x,p) 


Following steps of Lemma 11.35 of [FKOGj (which relies on Assumption 12.3p we get that 

inf Ho{x,p : e,C) < Ho{x,p). 

0<e<i4eC^{R) 

From the proof of Lemma B.IO in |FK06] . we have 


sup Hi{x,p:9,^)> inf liminft ^InE^ 
o<e<i,5eC2(R) 


Jo V{Ys’^;x,p)ds 


Thus, we need to show that, irrespective of the initial distribution, 


lim inf t ^ In Fi 

t^OO 




> Ho{x,p). 


The proof of this claim depends on the Assumption 12.31 We define the occupation measures of 
the process: 

Recall that P(M) is a separable metric space under the Prokhorov metric and that weak con¬ 
vergence of measures is equivalent to convergence in the Prokhorov metric. Let Qt,yo denote 
the probability measure on P(M) induced by the occupation measure pt of Y when Yq = y^. 
In other words, for A E .S(P(M)) (the borel sigma-algebra on P(M)), 


Qt^y^{A) = P{pt{-)€A\Yo = yo). 


Lemma 8. inf^gpj-R) lim inff_>.oo f ^ In 


AoV(Yf'^-,x,p)ds 


> Ho{x,p). 


Proof. Define (j) : V(R) ^ M by (/)(/r) = J V{y, x,p)p{dy). Take iq E P(K), and let R(Pi,r) 
denote the open ball in P(M) of radius r, centered at zq. Fix E 'P(M), then there exists a 
compact set K in M such that I'liK) > 0. The key ingredient in the proof is the uniform LDP 
lower bound for the occupation measures: 


lim inf — log 

t—>-oo t 


inf 

yo&K 


Qt,yo{B{vi,r)) 


> -J(Pi;x,p). 


(32) 


This is obtained from Theorem 5.5 in [DV83| under Assumption 12.31 While the statement 
of Theorem 5.5 in [DV83| is in terms of a process level LDP, by the contraction principle it 
ensures the uniform LDP lower bound (I32p for the occupation measures 
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We now compute 


lim inf - log 
t—>co t 


JlV(Y^’^-,x,p)ds 


= lim inf — log 
i—>-00 t 

- i l{yoeK} 

> lim inf — log 
t—>-oo t 


inf Ey^ 

yo&K \ J 


= lim inf — log 
i—>-00 t 


inf 


+ lim inf — log vi [K) 

t —^CXD ^ 


yoeK J^6P(r) 




> lim inf — log 
t^oo t 


1 


inf 

yo&K JfieB(i7i,r) 

> inf 0(/u) + lim inf — log 

lieB{i7i,r) t^oo t 

> inf 0(^) - J(z7i;x,p) 

tieB{ui,r) 




inf 

yo&K 


by (f32]l . By Lemma [9] (see Appendix), 4> is a lower semi-continuous function, and so < 

lim,._^o inf^gB{t 7 i,r) Thus taking limit as r —>• 0 we get 


lim inf — log E^ 
t^oo t 


Joy{yE^\x,p)ds 


> (/)(zA) - J(i7i) 


(note that since V is bounded below, > —oo, and so (()(ti) — J{ui;x,p) is well dehned 
and not —oo + oo). Since vi is arbitrary, we get 


lim inf — log Ei 

t —^OO ^ 


=/o V{Ys'’^\x,p)ds 


> sup {(^(hi) - J(z7i;x,p)}. 
iuG7^(R) 


This holds for every vi E P(M) and so 


inf lim inf — log E\ 
i^ieP{iR) t^oo t 


=/o V{YE^\x,p)ds 


> sup {(l){vi) - J{vi]X,p)]. 
ui eP(K) 


This concludes the proof of the Operator Inequality (j3ip . □ 

Remark 3.2. In the special case Y^’P also has a reversible invariant measure 7rP{x, •) we could 
also follow the arguments for Lemma 5-4 in jFFK10f using the Dirichlet form representation 
ofj (HHI). 

Proof of Lemma [T1 By Lemma [6] and Operator Inequality ()3ip . it follows that vf" is a sub¬ 
solution and vf- a supersolution of the Cauchy problem (I2TI1 : dtu{t,x) = HQ{x,dxu{t,x)) 
with u(0, x) = h{x). If the comparison principle holds for the Cauchy problem (I21h . then 
Lemma [3 gives us vf’ = and that ^ Uq = vf' = uniformly over compact subsets of 
[0, T] X M X R. □ 
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3.3 Exponential tightness 


Proof of Lemma [21 We prove exponential tightness using the convergence of and appeal¬ 
ing to supermartingale arguments (see Section 4.5 of |FK06] ). 

Let f{x) := ln(l -|- x^), so f{x) —)• oo as |x| —)• oo, and also ||/'|| -|- ||/''|| < oo. Define 
fe{x,y) := f{x) -|- eC(y) where Q is the positive Lyapunov function satisfying Assumption 12.41 
(with 6 = 1). Then, for any c > 0, there exists a compact Ac C M such that fe{x,y) > c, 
Vy G M, Vx ^ Kc- 

Observe that by (l2^ (with 6 = 1) 


HJ^{x,y) = ee 


< b{x,y)d^f{x) + -(T^(x,y)(9j,/(x))^ e bo{x,y)d^f{x) + -a‘^{x,y)dl^f{x 


+ {e 


f{x+ek{x,y,z),y) — f(x,y) 


— 1 — pk{x, y, z))ui{z)dz -f e 


By choice of /, growth conditions on the coefficients and Assumption 12.41 we get there exists 
C > 0 such that 

sup iLc/c(x,7/) <C< oo, Ve > 0. 

xGR^yGR 


Since e(/dW,t,T.t)-/dW,o,n,o))A-/o 
tional stopping 


is a non-negative local martingale, by op- 


P(Ac,t i 

7c(W,t,n,t) fe{xo,yo) 


< E 


exp 


- / HJ,{X,^s,Y,,s)ds 


< 1 . 


Therefore for each c > 0 


eIn P(W,t ^ Ac) <tC - /c(xo, yo) - c 

As C is fixed and independent of c (which we can choose), {Ac^ijoo is exponentially tight. 

Remark 3.3. A similar argument can be used to verify the exponential compact containment 
condition in Corollary 4-17 in \FKO0f . which would give us {Acc}e>o is exponentially tight. 

□ 


4 Examples 

4.1 Model for stock price with stochastic volatility 

We consider the stochastic volatility model for stock price suggested by Barndorff-Nielson and 
Shephard |B-NS0l] . Let Xt denote the logarithm of stock price and Yt the stochastic volatility. 

dXt = {r- ^Yt)dt + y^tdWt 

dYt = -^dt + dz}^\ 
d 
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where Wt is a standard Brownian motion and is an independent non-Gaussian Levy 

process with intensity jv{dz)\ the parameter 0 < <5 <C 1 denotes the mean-reversion time scale 
in stochastic volatility. The process Z is often referred to as the background driving Levy process 
(BDLP). If we are interested in pricing options on the stock which are close to matnrity, we will 
only be interested in small-time asymptotics of the model. We thus scale time by a parameter 
0 < e <C 1, to get 

dX,,t = e{r - + yTeV^tdWt 

2 ( 33 ) 

dY,^t = -^Y,,tdt + dzli\ 

The multi scale structure comes from the fast mean reversion in stochastic volatility and the 
small time to maturity. We are interested in the situation where time to maturity (e) is small, 
but large compared to mean-reversion time (5) of stochastic volatility. The interesting regime 
as seen in |FFK12j is when 6 = The generator of is given by: 

Yef{x,y) =€ (^{r - ^y)dj;f{x,y) + ^yd^J{x,y)^ 

+ ^ (^-ydyfix,y) + j {f{x,y + z)-f{x,y))u{dz)^, 

for / G C2(]R2). 

For this example, since the coefficients are x-independent, the perturbed operator is 
the same as Ci, the generator of Y: 

Yifiy) = -yf'iy) + J (/(y + z) - f{y)) u{dz), for / g C^{R). 


We can obtain the limiting Hamiltonian Flo by solving the eigenvalue problem (1161) . Here 
V{y,x,p) = V{y,p) = hyp^. Hq{p) is the eigenvalue A of the eigenvalue problem 


-yf{y) + J {f{y + z)- f{y)) y{dz) + -yp^f{y) = Xf{y). 

Note that f{y) = e 2 y and X{p) = J I e 2 ^ — 1 I u{dz) satisfy the eigenvalue problem. So 


Ho{p) = X{p) = f ye 2 ^ — ij i/[dz). In this example, in the absence of a Lyapunov function 

C satisfying Assumption 12.41 we give a slightly altered proof as follows. The following proof 
assumes Hq{p) is hnite. 

To verify Condition 13.11 for / G and 0 < 0 < 1, we define /o,e := f{x) + e{{l — 
0)g{x,y)+9({y)), where g{x,y) := ^{f'{x))‘^y (logarithm of the eigenfunction) andC(y) := C^y, 
C := supj, \f'{x)\. Then we get 

HJo^^(x,y) < Jy(/'(x))^ e((r - |)/'(x) -h Jy/"(x)) (1 - e)e~^c{ + Qe~^ll{ 


= (l-0)A(/'(x)) + 0 


-V[C^ - -Sr{x)f) + / (e^ ^ - \)v{dz) 


+ e((r- |)/'(x) + ^y/"(x)). 
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Thus satisfies Condition 13. fill Condition 13.1121 is immediate and 


lim sup e < inf 

e - j .0 ’ 0 < 6»<1 

< lim sup 
61-i.O 


(1 - 0)A(/'(x)) + 0sup -y{C^ --inx)f) + / - l^dz) 


(1 - e)X{f'{x)) + 0sup (-2/(C2 - ^{f'{x)f) + I- iMdz) 


= A(/'(x)) =: Hoif'ix)). 


Similarly, to verify Condition 13.21 define ■= fix) + e((l + 9)g{x,y) — 0({y)). It is 
unnecessary to verify any operator inequality as the limiting operators Hq and Hi coincide and 
equal Hq. 

Remark 4.1. Recall the definition ofV at the beginning of section \2. 1[ := Viy, x,p) + 

\boix, y)|+(T^(x, y). In general, in case we have a solution to the eigenvalue problem defining the 
Hamiltonian Hq, then the exact same proof as above using /o,e = fix)+ei{l — 9)gix,y)+9(iy)), 
with gix,y) the logarithm of the eigenfunction and C satisfying the V-multiplicative ergodicity 
condition 

e~^£^’^e^iy) < —cV^{x, y) + d, for c> l,d > 0 

is enough to conclude our large deviation results (provided V has compact finite level sets, as 
it was above). 


In Barndorff-Nielsen and Shephard |B-NS01| . the BDLP, Z, is assumed to have only positive 
increments. A simple example of such a Levy process is a jump process taking finitely many 
jumps that is the Levy measure is ^{zi) > 0 where Zi > 0, i = 1,2,..., k. We can then explicitly 
compute Hq(p) and its Legendre transform L(p). As seen in [FFK12] (Lemma D.l in [FFK12] ). 
since Hq{p) is not state dependent, we get the rate function to be I{x,XQ,t) = tL ()• In 
finance, a common example is where Z is a gamma process, in which case v{dz) = ^e~^^dz, 
a,b > 0. Then 


Hq{p) 


a In (^1 + if - <p < 

oo if > 2b, 


and the rate function is given by I{x‘,XQ,t) = tL ( ^°^ , where 


l{q) 


—a + + 2bq‘^ — a In 26 + a In + 2bq^^ 

< 0 

—a — a? + 2bq^ — a In 26 + a In - + 2bq‘^'j 


if g > 0 
if g = 0 
if g < 0. 


This rate function then gives the asymptotic behavior of a European Call option on the stock. 
Let K denote the strike price and then for Sq = < K (out-of-the-money call), 


lirn e log E - K]^ = -/(log A; xq, t), 


where maturity time T = et. This follows from Corollary 1.3 in [FFFlOj . 
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4.2 Model for self-regulating protein prodnction 

The simplest model for translation of protein from DNA is the system below, with a gene that 
is either in its “on” state Gi, or in its “off’ state Go, and in which the protein activates the 


changes from “ofF’ to 

“on” state: 






(1) 

Go + P ^ 

Gi+P 

(3) 

Gi 


Gi + P 

(2) 

k' , 

Gi + P 

Go + P 

(4) 

P 

_1 

0 


Suppose the amount of protein P is of order 1/e, whose rate of production ^2 = l/eK 2 , 
while its rate of degradation k'^ = K 3 ; where K 2 , are of 0(1). The amount of genes in the 
“on”- and “ofF’-state is G {0,1}, their total amount always equaling 1, and suppose the rates 
of changes of the gene from the “on”-state to the “off’-state and back are very rapid due to 
its regulation by the large amounts of protein k[ = ki,k'_^ = k_i, where ki,k_i are of 0(1). 
This system is characteristic of eukaryotes, where the gene switching noise dominates over the 
transcriptional and translational noise. We can represent the changes in the system using the 
process Xf, for the count of protein molecules normalized by e, and 1 / for the (unnormalized) 
count of “on”-gene molecules. A diffusion process is a good approximation for the evolution 
of as long as the count of proteins is not too small, that is, the unnormalized count is 
e and X^ ~ 0(1) f |KKP12] gives a rigorous justification of diffusion approximations for 
Markov chain models that apply in stochastic reaction kinetics). This diffusion solves dX^^t = 
b{X^^t,Ye,t)dt -|- y/ea{Xf:^t,Y^^t)dWt with drift b{x,y) = K 2 y — (protein production has only 
two possible values: it will be 0 when y = 0 , or K 2 when y = 1), with diffusion coefficient 
a‘^{x,y) = K 2 y + K 3 X, and initial value X^^fl = xq > 0. Changes in the amount of proteins due 
to other independent sources of noise, such as errors after cell splitting, can be modelled by 
an additional jump term for where the jump measure i'i{dx) can be as simple as vi{z) = 
+ ^ 6 +i{z), producing 

dX^^t = iK2Ye,t - H3Xe,t)dt + 6(^2d/,* -F K 3 X^^t)dWt + € J lx,>eZN'^ {dz, dt) 

The amount of genes Gi in the “on”-state is a rapidly fluctuating two-state Markov chain 
Y on {0,1} with rates ro_).i(x) = \kix and ri^Q{x) = \k-\x that depend on the normalized 
amount of protein (note that the amount of genes Go in the “ofF’-state is 1 — Y). This chain 
is reversible, and for each x > 0 it has a unique stationary distribution 7 r^(l) = 1 — 7 r*( 0 ) = 
ki/{ki + K_i). 

Signalling proteins such as morphogens have to be in the right range of concentrations to 
avoid triggering the expression of genes at the wrong times. The probabilities of their amounts 
being out of range are given by the Large Deviation Principle for as e —)• 0, for which we 
need to obtain the solution to the eigenvalue problem for the operator V{y;x,p) + CP where 

CYf{y) = ro^i{x){f{y + 1) - f{y))ly=o + ri^o{x){f{y - 1 ) - f {y))ly=i. 
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In order to solve {V{y,x,p) + = Ae“i for A, let = ai(x), = ao(x), for 

some 01,00 strictly positive functions. Then 

(k 2 - K3x)pai(x) + (k 2 + K3x)p^ai(x) + + e~^ - 2 )oi(x) + k_ix(oo(x) - ai(x)) = Aai(x) 

-K3xpao(x) + K3xp^ao(x) + + e~^ - 2 )ao(x) + kix(oi(x) - ao(x)) = Aao(x) 


equivalently, with o(x) = oi(x)/oo(x). 


1 


(k2 - K3x)p + (k2 + K3x)p^ + K-ixf—— - 1) = -K3XP + K 3 Xp^ + Kix(o(x) - 1) 

o(x) 

which, since o(x) has to be positive, gives 

, , -B + \/52 _ 4AC , 2 . N ^ 

o(x) =-—-, yl = Kix, B = —K 2 P — K 2 P + (k_i — Kijx, C = —K_ix 

and consequently, using notation above, 

Ho{x,p) = -K 3 XP + K 3 xp^ + Kix(a(x) - 1) + - 2) 


Note that when k_i = ni then 


K2p{l + p) + V {k2p{ 1 + p)f + (2kix)2 

ox =-- 

2kix 

and 

Ho{x,p) = -K 3 p {1 - p)x + ^K 2 p {1 + p) + ^\/ {K 2 p{l+p)f + (2kix) 2 - Kix + + e“^ - 2). 

Note that Hq{x,p) is both convex in p and continuous in x. 

If one were to use an approximation of the evolution of the normalized protein amount 
by a piecewise deterministic process then (without additional noise) 

^^PDMP ^ 

while is the same fast Markov chain on {0,1}. In this case V{y;x,p) = (k 2 — H 3 x)p and the 
Hamiltonian (when ki = k_i) becomes 

Hq^^^^{x,p) = -K3PX + ]^K2P+ ]^\/{f^2PY + (2kix)2 - KiX. 

which is easy to compare to the Hamiltonian Hq of the diffusion process taking into account 
the small perturbative noise arising from randomness in the timing of chemical reactions and 
from randomness in the outcomes of cell splitting. 
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A Appendix 


Lemma 9. Fix G M and let (j) '■ 'P(IR) —)• K 6e defined by (pifi) = J V{y, x,p)iJ,{dy). Then, 
(j) is a lower semi-continuous (l.s.c.) function onP(M). 

Proof. For the rest of the proof, we will write V{y) for V{y]x,p). Let Vm ■= V ■ lv<M + M ■ 
for M > iniyV{y). To show that 4>{iJ,) is l.s.c, it is sufficient to show that if —)■ p 
weakly, then (p{p) < liminfn^oo 0(/Un)- Assume pn —^ P weakly. Then 



lim 

n—^oo 



by definition of weak convergence of measures, since Vm is a bounded function. By the mono¬ 
tone convergence theorem we get 


4>{t) = Vdp = lim / VMdp 
I M—>-oo / 


= lim lim / VMdpr, 

M—^oon^oo J 

= sup lim / VMdpn 
M n^ooj 


< lim inf sup / VMdpr, 
M J 


= lim inf / Vdp„ 
n^oo I 


by Monotone convergence theorem 


= lim inf 4>{pn) 

n—^oo 


□ 


Lemma 10. Let ui be a bounded, upper semicontinuous (u.s.c.), viscosity subsolution and U 2 
a bounded, lower semicontinuous (l.s.c.), viscosity supersolution of dtu{t,x) = HQ{x,dxu{t,x)) 
respectively. //ui(0, •) < rt 2 ( 0 , •), and Hq is uniformly eontinuous on compact sets, then ui < U 2 
on [0,T] X M for any T > 0. 


Proof. Suppose 

sup {ui{t,x) — U 2 {t,x)} > A > (5 > 0. 

t<T^x 

Let g{t, x) = ln(l -|- x^) -j- t^. Define 


(A.34) 


fi{t,x,s,y) = ui{t,x) 


U2{s,y) - 2 



fi {g{t,x) +g{s,y)) - At. 


Fix /3 > 0 and let {fi, x^, s^, yfi) denote the point of maximum of V’ in ([0, T] x M x [0, T] x R) for 
e > 0. Since ui,U 2 are bounded, for fixed /3 > 0, there exists an > 0 such that \xe\-, \ye\ < Rp 
for all e > 0. 
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Using 


we get 
1 


In 1 + 


•0(4, Xe,te,Xe) + Ve, S^, Ve) < 20(4, y^), 

\Xe - Ve? + |4 - 


< A(Se - 4) + Ui{t^,Xe) - Ui{Se,ye) + U2{te,Xe) - U2{Se,ye) 

< 2AT + 2||ui|| + 2|tt2|| =: C < oo. 


which gives us 

\xe - yep + |4 - 4p < ee^*^- 

Therefore |xe — ye|, |4 — 4| —>■ 0 as e —>■ 0. 

Let 


0 i(t,x) := U2(se,ye) + ^ In ( 1 + 


X - ye\ + 0 - S, 


5 |2 


and 


02 (s,?/) := Ui(4,Xe) - - In ( 1 + 


ke - 2/P + |4 - sl 


+ H9{t,x) + g{Se,ye)) + At 


- /? i9ite,Xe) + y(s, y)) - Ate. 


Then (4,Se) is a point of maximum of ui{t,x) — 0i(t,x) and {se,ye) is a point of minimum of 
U 2 {s,y) — 02(s,y). Since ui and U 2 are sub and super solutions respectively, by the definition 
of sub and super solutions we get 


U-Se 




1 + 




and 


te-Se 


1 + 


he-^eP + Pe-^eP 


~\~ A-\- 2 fjte ^ Hq [ Xe, 


- 2l3Se > -ffo 2/e 


+ 


2/3Xe 


^ |x-e-yeP + Pe-^d^ 1 + x^ y ’ 

2/3ye \ 


Xe-Ve 


Subtracting (|A.36I) from (IA.35p . we get 
A+2/3(4+4) < Hq i Xe, 


+ 


]_ _|_ |Xe-yeP + Pe-5ep 1 + y- 

2/3x 


r,2 


Xe-Ve 


I he-yeP + Pe-^eP 1+X 


^2 


-Ho ye, 


(A.35) 

(A.36) 
2^ye ^ 


^ , |Xe-yeP + Pe-^eP 1 + y^ j ' 


(A.37) 

Since i7o(',') is uniformly continuous over compact sets, and since |xe — ye| —)• 0 as e —>• 0 
(for fixed /3), the right-hand side of the above inequality goes to 0 as e —>• 0 and /3 —>• 0 (note 

—ye _ „ _ 

that the terms —— _ ,%, ^ and are bounded and that IxeLi/el < Rb for each 

, \xf-yf\.‘+\tf-Sf\.‘ ’ 1+x; 1 + 1 /; I tijwei _ p 

/ 3 ). 


1 + 


Taking e —)• 0 and then /3 —)• 0, we get 


A < 0, 

which contradicts ()A.34p . Therefore we must have 

SUp{Ml(t, x) — U 2 {t, x)} < 0 

t,X 

which gives us ui < U2. 


□ 
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